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Abstract 


Direct  solid-liquid  interfacial  energies  are  difficult 
to  measure  experimentally.  A  relationship  between  the  solid- 
liquid  surface  tensions  of  the  basal  and  prismatic  planes 
of  ice  at  -38.8°C  is  calculated  using  equilibrium  thermodynamics 
and  homogeneous  nucleation  theory  for  a  nucleus  shape  based 
on  that  of  a  unit  cell.  Using  contact  angle  theory  and 
experimental  data  the  solid-liquid  surface  energy  for  the 
basal  plane  is  calculated  as  17.0  dyne/cm  and  that  of  the 
prismatic  plane  as  28.9  dyne/cm. 
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1.0  Introduction 


The  surface  energy  of  an  interface  is  the  work 
required  to  create  a  unit  area  of  the  interface.  The 
following  four  surface  energies  can  be  defined  for  any 
liquid-solid-vapour  system. 

LV 

1.  -  the  liquid-vapour  interfacial  tension, 

2.  Y  -  the  solid-vapour  interfacial  tension, 

3.  Y  -  the  solid-liquid  interfacial  tension, 

Sb 

4  .  Y  -  the  grain  boundary  energy. 

In  this  work  the  system  considered  is  an  ice-water-water 

.SL 

vapour  system  and  values  for  5  the  solid-liquid  inter¬ 
facial  tensions  are  calculated. 

The  solid-liquid  interfacial  tension  is  important 
to  many  processes.  The  temperatures  at  which  solids 
nucleate,  whether  it  be  on  foreign  particles  as  in  heter¬ 
ogeneous  nucleation,  or  by  the  spontaneous  growth  of 
water  molecules  combining  together  to  form  an  ice  nucleus, 
homogeneous  nucleation,  are  largely  determined  by  the 
liquid-solid  interfacial  tension.  The  morphology  of 
crystal  growth,  particle  engulfment,  frost  heave  and 
pipeline  freezing  are  all  examples  of  phenomena  where 
surface  tension  has  a  critical  role. 

Of  the  previously  defined  surface  energies  only  the 
liquid-vapour  surface  energy  is  relatively  easy  to 
measure.  The  Wilhelmy  plate  method  is  one  well-documented 
method  (Davies  and  Rideal,  1963).  Essentially  the 
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method  measures  the  additional  pull  on  a  thin  plate  when 
it  becomes  partially  submerged  (see  Figure  1.1).  The  addi¬ 
tional  force  is  equal  to  the  product  of  the  perimeter 
of  the  plate  and  the  surface  tension.  Other  common  methods 
include  the  small  ring  method,  the  sessile  drop  method, 
the  drop  weight  method  and  pendant  drop  method, all  of 
which  are  documented  by  Matijevic  (  196  9  :  101-152).  The 
liquid-vapour  interfacial  tension  for  water-water  vapour 
is  calculated  in  chapter  three  of  this  work. 

The  ice-water  vapour  interfacial  tension  is  more 
difficult  to  measure  and  estimates  for  this  value  vary. 
Mason  (1952)  and  McDonald  (1953)  calculated  the  number 
of  hydrogen  bonds  that  must  be  broken  to  create  an  ice  surface. 

CD  £ 

At  0  C  their  estimates  for  range  from  85-112  dyne/cm. 

Ketchum  and  Hobbs  (1969)  used  a  direct  experimental  method 

involving  measurement  of  dihedral  contact  angles  and  cal- 

+  sv 

culated  a  value  of  109-3  dyne/cm.  These  values  for  ^ 

must  be  considered  average  values  since  they  are  not  spec¬ 
ific  to  a  particular  plane  of  ice. 

SL 

Methods  for  measuring  °  ,  the  solid-liquid  surface 

tension  are  numerous  and  varied.  Homogeneous  nucleation 
theory  (explained  in  more  detail  in  chapter  three)  and 
experiments  are  used  by  many  researchers  to  estimate  this 
value.  Turnbull  (1950)  was  one  of  the  first  to  use  homo¬ 
geneous  nucleation  theory  to  calculate  the  interfacial 
energy  of  solid  nuclei  and  their  corresponding  liquids. 
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Figure  1.1  Wilhelmy  Plate  Method 


The  pull  on  a  verticle  plate  per 
is  V  *cos  ©.  (Adapted  from  Davies 
1963:46) 


unit  perimeter 
and  Rideal, 
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A  summary  of  Turnbull's  and  other  experimenters'  estimates 

,SL 

for  2f  is  given  in  Table  1.1.  Other  methods  include 
theoretical  estimates,  Table  1.2,  direct  measurements 
of  grain  boundary  groove,'*'  measurement  of  dihedral  and 
contact  angles,  Skapski's  capillary  cone  method  (Table 
1.3)  and  crystal  growth  experiments  (Table  1.4).  As  seen 

,SL 

from  Tables  1.1  to  1.4,  the  values  for  o  for  ice-water 
vary  from  6  to  50  dyne/cm.  Note  that  only  Hillig  (1958) 
and  Simpson  et  al  (1974  a  and  1974  b  )  report  any  orienta¬ 
tional  dependence  for  the  surface  tension.  Although  Hesstvedt 
(1967)  realized  that  each  different  kind  of  face  on  a 
solid  nucleus  should  have  a  unique  value  for  surface  ten¬ 
sion,  he  did  not  calculate  them.  He  incorporates  a  geomet¬ 
ric  factor  which  varies  for  different  crystal  shapes 
into  his  equations.  The  value  he  obtains  is  thus  a  weighted 
mean.  The  approach  adapted  in  this  work  is  definitive 
because  it  maintains  the  identity  of  the  surface  tensions 
corresponding  to  different  faces  of  the  solid  nucleus. 

Since  "ice  generally  crystallizes  in  hexagonal  prisms, 
this  form  is  usually  taken  to  be  the  form  of  the  germ" 
(nucleus)  ( Duf our-Def ay ,  1963:162).  In  the  following 
chapters,  a  nucleus  shape  based  on  the  'unit  cell'  for 
ice  is  used,  ie  a  hexagonal  prism.  It  is  important  to 

■*"When  a  crystalline  solid  is  exposed  to  a  vapour  or  liquid 
grain  boundary  grooves  form  at  the  interface.  Measurement 
of  the  grain  boundary  groove  angle,  the  use  of  Young's 
equation  (see  chapter  three)  and  experimental  data  for  # 
permits  the  calculation  of  lfi/s  and  H  (see  Hobbs,  1974) 
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Table  1.1  Estimates  of  0  Based  on 
Nucleation  Experiments 


( Adapted 

from  Kallungal  (1975:62) 

Reference 

dyne/cm 

Wood  &  Walton  (1970) 

27.53  -  .  037 

Turnbull  (1950) 

32 . 1 

Carte  (1959) 

16 . 0 

Jacobi  (1955) 

15.5 

Hesstvedt  (1964) 

31.7 

Table  1.2 

S  L 

Estimates  of  H  Based  on 

Theories  (Kallungal  (1975:61) 


Reference 

dyne/ cm 

Defour  &  Defay  (1963) 

10-49 

Fletcher  (1959) 

20 

Turnbull  (1965) 

38-50 

Skapski  (1956) 

27 . 8 
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Table  1.3  Estimates  of  ^  Based  on 
Estimated  Measurements 
(Adapted  from  Kallungal  (1978:62) 


Reference  dyne/cm 


Jones  ( 1973  ) 

44 

-  10 

Jones  and  Chadwick  (1971) 

41 

±9 

Ketchum  &  Hobbs  (1969) 

33 

±  3 

Suzuki  &  Kuroiwa  (1975  ) 

45 

-  15 

Jones  (1970) 

46 

Skapski  (1957) 

44 

o 

i — i 

+  i 

s  u 

Table  1.4  Estimates  of  Based  on 

Crystal  Growth  Experiments 
(Adapted  from  Kallungal  1975:63) 


Reference 


Fernandez  (1967  a) 
Poisot  (  1968  ) 
Vlahakis  ( 1972  ) 
Hardy,  Coriell  and 
Hardy  and  Coriell  ( 
Kotler  and  Taishis 
Simpson  et  al  (1974 
Simpson  et  al  (1974 
Hillig  (1958) 


31 

49 

33 

Sekerka  ( 1971) 

25 

1968  ) 

16 

( 1968 ) 

20 

a) 

28 

b) 

6  . 

6  . 

iyne/cm 

.8  -  1.8 

.  3 

.2-4 

.0 

.0 

±  2 

.0  (prismatic  plane) 
24  (basal  plane) 

4  (basal  plane) 
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note  that  only  the  shape  of  the  unit  cell  is  chosen,  not 
the  relative  size.  For  example  the  nucleus  could  be  flat 
and  disk-like  or  long  and  columnar.  Equilibrium  theory 
is  used  to  define  each  surface  tension  in  terms  of  the 
measurable  quantities  temperature  and  pressure.  Stability 
analysis  will  be  used  to  illustrate  the  concept  of  a 
"critical  size"  nucleus  for  a  given  temperature  and  pres¬ 
sure.  Homogeneous  nucleation  rate  theory  and  experimental 
data  are  used  to  generate  a  relationship  between  the  sur¬ 
face  tensions  for  the  different  planes  of  the  nucleus. 

To  calculate  the  individual  values  of  the  surface  tensions, 
recent  work  on  contact  angle  theory  and  experiments  are 
used  to  calculate  the  solid-liquid  surface  tension  of  the 
basal  crystallographic  planes  of  ice. 


2 . 0  Theory 


This  chapter  develops  the  theory  necessary  to  arrive 
at  an  expression  for  the  equilibrium  size  of  a  solid 
nucleus  surrounded  by  a  liquid  in  terms  of  the  measurable 
quantities  of  temperature,  T,  and  pressure,  P.  The  postul- 
atory  formulation  of  reversible  thermodynamics,  as  treated 
by  Callen  (1960)  is  utilized  to  analyze  the  thermodynamic 
equilibrium  and  a  solid  nucleus  in  a  surrounding  liquid. 

The  chapter  is  divided  into  the  following  3  subsections: 

2.1)  Equilibrium  Conditions. 

2.2)  Stability  Theory. 

2.3)  Expression  for  Equilibrium  Size  of  Nucleus. 


8 


9 


2.1  Equilibrium  Conditions 

In  this  section  the  conditions  under  which  a  solid 
nucleus  is  in  thermodynamic  equilibrium  with  its  liquid 
phase  are  developed.  The  agreement  of  predictions  of  equil¬ 
ibrium  thermodynamics  with  our  intuitive  conceptualizations  of 
temperature  and  pressure  are  often  illustrated  (Callen,  1960). 
Unfortunately  this  cannot  be  said  for  surface  tension.  There¬ 
fore  the  opportunity  to  illustrate  that  the  behaviour  of  the 
partial  differential  of  internal  energy  with  respect  to  area 
agrees  with  our  intuitive  notions  of  surface  tension  is  taken. 

To  illustrate  the  agreement  of  equilibrium  thermodynamics 
to  our  intuitive  concept  of  surface  tension,  the  system  shown 
in  Figure  2.1  is  considered.  Single  component  films  A  and  B 
each  float  on  an  inert  liquid.  They  are  separated  by  a  movable 
inert  barrier.  The  system  is  in  equilibrium  at  constant, 
uniform  temperature.  The  energy  minimum  principle  states 
that  dUj-,  the  total  internal  energy  change  for  the  system  is 
zero : 

JUT  -  O. 

(2.1.1) 

S’ 

The  internal  energy,  U  ,  for  the  surface  layers  A  and  B 

s  s 

is  a  function  of  the  entropy,  S  ,  the  area,  A  ,  and  the  number 

£ 

of  moles  of  material,  N  ,  in  the  surface  layer: 


u3--  us(s*  AS,NS) 


(2.1.2) 


Figure  2.1  Single-Component  Films  A  and  B. 
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When  any  two  phases  are  in  contact  with  each  other 
they  are  separated  by  an  interfacial  layer  see  Figure  2.2. 
The  value  of  a  given  property  will  vary  continuously  from 
its  value  in  the  liquid  phase  to  that  in  the  solid  phase. 

In  the  system  to  be  discussed  in  the  following  pages  the 
surface  layer  is  replaced  by  a  dividing  surface. 

At  this  surface  there  is  a  discontinuity  in  the  value 
of  the  property.  This  is  a  conceptualization  conceived 
by  Gibbs  (1961)  while  studying  the  phenomena  of  capillarity. 
The  position  of  this  dividing  surface  is  defined  as  being 
where  the  adsorbtion  of  component  one,  "the  solvent",  is 
zero  where 


5 


(2.1.3) 


are  commonly  known  as 


'excess  properties'  in  that  they  represent  the  difference 
between  the  total  amount  of  property  and  the  amounts  in 
the  liquid  (')  and  solid  (")  phases. 

They  are  commonly  defined  as  follows: 


t/=  U  -  u'~  u 


(2.1.4) 


S  -  s '  -  s " 


(2.1.5) 


I  I 


and 


(2.1.6) 
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Figure  2.2  Variation  of  a  Property  in  the  Surface  layer. 

- actual  variation 

- modelled  variation 
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Applying  equation  2.1.2  to  the  system  in  Figure  2.1  yields 


du‘ *</<//  =  auras'  +  -h  djxA  J asa 

&  s„  At  s  a  )ct  u* 

3<4) 

^.<^3 


^4 

.S 


5 


*h 


(2.1.7) 


+  3  U.  ]  d  f  ^Jaa(b\cI^s  ' 

3’  "<*  ^  j s* 1  d  ^0  As Q  >  M 3 

In  equation  2.1.7  the  partial  differentials  are  more 
commonly  known  as  the  following  intensive  parameters: 


A,  /V 


temperature,  (2.1.8) 


au )  = 

and  ^  ^  ^  '  the  chemical  potential.  (2.1.9) 

The  following  discussion  illustrates  how  the  partial 
differential,  ( c)  U  /d^)s  ^  ,  behaves  in  a  manner  associated  with 
'$  ,  the  surface  tension. 

The  following  constraints  are  evident: 


1) 


dU„  =  O) 


2) 


d  N8-  ° 


3) 


„  s  5 

^R  h  ^  8 


=  constant 


/ 


implying 


d  s„  ^  i 


B  > 


4) 


f\l  J 


g  =  constant 


implying 


J  /)  -  -  d  $ 


e 


(2.1.10) 

(2.1.11) 

(2.1.12) 

(2.1.13) 

(2.1.14) 

(2.1.15) 
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Using  equations  2.1.1  and  2.1.8  to  2.1.15  to  simplify 
equation  2.1.7  yields: 


(\  -Ja)ois, 


d 


A«  '  d  J 


(2.1.16) 


5  S 

The  variables  dSA  and  dAA  must  be  permitted  to  vary 
independently  of  each  other.  This  leads  to  the  two  equilib¬ 
rium  conditions: 


8 


and 


-  c)  ii 


B 


Mi  ' s’ ,  n!  spia/saJNs 


(2.1.17) 


(2.1.18) 


P  >  Nft 

Physically,  at  equilibrium,  the  surface  tensions,  2f  ^ 


and  ,  are  expected  to  be  equal;  just  as  the  quantities 

in  equation  2.1.18  are  equal. 

To  examine  the  effect  of  a  perturbation  on  the  system, 
consider  the  effect  of  imposing  a  greater  surface  tension 
of  film  A  than  film  B  on  the  system  illustrated  Figure  2.1. 
This  state  of  the  system  will  be  referred  to  as  state  1  and 
the  equilibrium  state  will  be  referred  to  as  state  2. 

Figure  2.3  is  an  illustration  of  state  1.  Note  that 
the  temperatures  of  films  A  and  B  are  equal,  only  the  surface 
tension  has  been  altered. 

The  total  energy  U,  for  each  state  is: 


(/)  0 )  O) 

U  --  Ufi  *  u 


B 


(2.1.19) 


( z >  60  Ci) 

U  --  Ua  I-U 


6 


and 


(2.1.20) 


. 


maaasassssi 
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e 


Film  B 


A 


Inert  Liquid 


Single-Component  Films  A  and  B, 

Y/> 


where  'f  is  greater  than  Yd1. 


Figure  2 . 3 


LiSBasS 
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The  change  in  energy,  A  U,  between  states  1  and  2  is 


(V)  //) 

A  U  =  U  -  U  .  (2.1.21) 

Substituting  equations  2.1.19  and  2.1.20  into  equation  2.1.21 
yields:  (i)  (,)  .  ,  (2.)  a  )  \ 

)  +  (  U8  -Us  )  . 

(2.1.22) 

A  Taylor  series  expansion  about  the  equilibrium  point, 
state  2,  exclusive  of  terms  higher  than  first  order  yields: 

C/i  cj) 


O)  ^  .  (,)«.  .  I/'  % 

Ur  :  U a  f  '  Aa  +-  • 

^ h) 

(CCh.. 


c/l 

^6  *  ^8  f 


c)  U 

c)  Ag  ) 


(2.1.23) 


(2.1.24) 


These  can  be  simplified  to: 


and 


{/•> 

if) 

V* 

^  , 

(2. 

.1.25) 

[  1 1 

i  Cl) 

y 8 

aA8  . 

(2. 

.1.26) 

and 

2.1.26  into 

equation  I 

> .1.22 

(i ) 

u  ) 

yields : 

A  U-  /a  A'!’  -y  aA6  -  (2.1.27, 

Recalling  constraint  2.1.14  and  substituting  into  equation 
2.1.27  yields  the  following  expression: 

( i)  i  ,)  \  .  ()  ) 


AU--  -  (  -  r3 


-JL')aa 


A  • 


(2.1.28) 


17 


The  energy  minimum  principle  implies  that  A  U  is  a  negative 


quantity . 


U)  M 


Since  is  imposed  as  being  a  positive  quantity, 


u) 


A  ,  must  be  a  positive  quantity.  This  is  in  agreement  with 
our  intuitive  feeling  that  the  film  A  with  a  greater  surface 
tension  increases  the  area  it  occupies  at  the  expense  of 
film  B.  Physically  it  is  expected  that  the  film  with  the 
larger  surface  tension  will  expand  its  area. 

In  a  laboratory  situation  intensive  parameters  are  often 
the  most  readily  varied.  It  is  simpler  to  measure  an 
intensive  parameter  like  temperature  than  it  is  to  measure 
S,  the  entropy.  The  partial  Legendre  transform  that 
replaces  temperature  for  entropy  as  the  independent  variable 
in  the  energy  representation  is  the  Helmholtz  free  energy, 

F.  The  Helmholtz  free  energy  is  a  function  of  temperature, 
volume  and  mole  fractions  N()^  ...A/p: 


F  r  h  (  l}V,  Nf,  A /*>••• 

The  complete  differential,  dF,  is, 

d F- - S  d  T  -  /°dVf  d  NL  l  r  »,  . . .  r  ( 2 .  i .  3  0 ) 

The  equilibrium  condition  in  the  Helmholtz  free  energy  repres¬ 
entation  is: 

"The  equilibrium  value  of  any  unconstrained  internal  parameter 
in  a  system  in  diathermal  contact  with  a  heat  reservoir 


Nr'* 


(2.1.29) 
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minimizes 
(equal  to 


the  Helmholtz  potential  at  constant  temperature 
that  of  the  heat  reservoir)"  (Callen  1960:105). 


Thus 

subject  to 


df  =  ot 

the  auxiliary  condition  that 


(2.1.31) 


T  -  T  r 

c 

where  T  is  the  reservoir  temperature.  The 
energy  is  the  representation  that  is  most  c 
the  forthcoming  discussion  of  a  solid  nucle 
with  a  liquid. 

Consider,  now,  the  thermodynamic  equil 
nucleus,  assumed  to  be  under  hydrostatic  st 
Figure  2.4,  surrounded  by  a  liquid,  each  phas 
facial  zone,  consisting  of  a  finite  number, 
ing  components.  This  system  is  isolated  by 
maintain  a  constant  temperature.  The  total 
is  held  constant,  is  the  sum  of  the  nucleus 


(2.1.32) 
Helmholtz  free 
onvenient  for 
us  in  equilibrium 

ibrium  of  a  solid 
ress,  as  shown  in 
e,  including  inter- 
r,  of  non-react- 
a  heat  bath  to 

volume,  V  ,  which 

,  il 

volume,  V  ,  and 


/ 


the  liquid  volume,  V  .  The  volume  of  liquid  phase  is  far 
greater  than  that  of  the  solid  phase  nucleus  making  changes 
of  pressure  due  to  growth  or  dissolution  of  the  solid  nucleus 
negligible.  As  mentioned  earlier,  under  these  constraints 
the  appropriate  potential  to  express  the  equilibrium  condition 
is  the  Helmholtz  free  energy,  equation  2.1.30.  The  condition 
of  equilibrium  for  this  system  is  that  the  total  change  of 
free  energy  must  equal  zero, 
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Figure  2.4  Thermodynamic  Equilibrium  of  a  Solid  Nucleus 
where 

(  ) denotes  liquid  phase 
(  11  )  denotes  solid  phase 
(s  ) denotes  interface. 
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(2.1.33) 

where  N  denotes  the  mole  numbers,  of  each  component,  N  j  , 
Nj  ,  ..  .  ,  N  c  •  This  equilibrium  condition  is  subject  to 
the  following  constraints: 

i)  the  total  volume  is  constant,  ie. 


(2.X. 34) 


ii)  the  total  number  of  moles  of  each  component  is  constant 


le 


constant 


(2.1.35) 


To  maintain  the  generality  of  the  problem  it  is  assumed 
that  the  solid  nucleus  has  c*  plane^  surfaced  sides.  Figure 
2.5  represents  the  growth  of  one  such  plane.  is  the  area 

of  this  particular  face,  h £  is  the  perpendicular  distance  from 


the  face  to  the  centre  of  the  nucleus 


V.. 


is  the  original 


n 


volume  of  the  solid  nucleus  and  dV  is  the  volume  of  the 


addition  to  face  A 


*  ' 


V  is  the  volume  of  the  liquid  and  dV 


(the  volume  of  liquid  solidified)  is  considered  to  be  negligible 
compared  to  the  total  volume  of  liquid. 

The  total  change  in  the  Helmholtz  free  energy  is  the 
sum  of  the  change  in  the  free  energy  of  the  bulk  phases  and 
the  change  in  free  energy  of  the  surface,  ie: 


dF=  dF'+dP'V  d  Fs  +  JF* 

_1This  paper  will  be  restricted  to  plane  surfaces. 


(2.1.36) 
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Figure  2.5  Growth  of  One  Face  of  a  Crystal  (Adapted  from 
Defay  and  Prigogine,  1966) 
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The  change  in  free  energy  of  the  fixed  volume  is  zero. 


Therefore 


d  P*  -  O  . 


(2.1.37) 


The  first  term  in  equation  2.1.36,  dF  ,  is  the  change  of 
free  energy  in  the  volume  V  '  ; 


d  F ^  -~F  cA  V  +  ^  d  A/ ^  7 


(2.1.38) 


whe 


th 


re  yCf i  is  the  chemical  potential  of  the  i  component 


/  Uh 

and  dN ^  is  the  change  in  the  number  of  moles  of  the  i 

/ 

component  in  the  volume  V 

The  change  in  free  energy  due  to  the  change  of  phase 
of  the  volume  dV  is 


r*  /  i 

1 1  #  «  I f  •  /  / 

dp  -“P  d  V  +  o/hfc  . 


(2.1.39) 


dF  ,  the  change  in  free  energy  for  ,  planar  surfaces  is 


d  F  =  ^  2T  d  -h  £  JA/- 

Using  the  relationship, 


s 

‘K  ' 


(2.1.40) 


dv'tdv  -  o. 


(2.1.41) 


substituting  equations  2.1.37  to  2.1.40  into  equation  2.1.36 
and  re-arranging  the  terms  yields: 


dF:  -  (  P"- p‘)c iv"f  +• 

- 


i  c/N-  +  if'  +  ^ 

1,-1  u  -  /  y 


(2.1.42) 
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The  energy  minimum  principle,  equation  2.1.33,  with  con¬ 
straint  equations  2.1.34  and  2.1.35  applied  to  equation  2.1.42 
gives  two  distinct  sets  of  equilibrium  conditions:  one 
involving  surface  tensions  and  the  other  chemical  potentials. 

Since  moles  and  volume  must  be  able  to  vary  independently 
it  is  possible  to  divide  equation  2.1.42  into  the  following 
two  equations: 


(2.1.43) 


The  following  discussion  will  first  treat  equation  2.1.43 
and  then  equation  2.1.44. 

From  Figure  2.5  the  following  geometric  relationship 
is  evident: 


d  V  -  Af  d  hf  • 


(2.1.45) 


Substituting  equation  2.1.45  into  equation  2.1.43  yields: 


(2.1.46) 


From  geometrical  considerations  it  can  be  shown  that; 


Substituting  equation  (2.1.47)  into  equation  (2.1.46) 


and  collecting  terms  gives: 


- 


7 


('  Xk  t-  1/)^ 

-  «  Jbp 


(2.1.48) 


where 


A?  (p"-p ') 


(2.1.49) 


For  a  particular  face  the  solution  to  this  equation  is 


(2.1.50) 


This  form  of  solution  is  known  as  the  Wulff  equation 
(Defay  and  Prigogine,  1966:296), 


(2.1.51) 


It  can  be  shown  (Defay  and  Prigogine,  1966:296)  that  "in 
every  crystal  which  is  in  equilibrium  there  exists  a  point 
0  such  that  the  Wulff  relations  are  satisfied.  The  Wulff 
form  is  therefore  the  only  equilibrium  form"  0,  is  defined 
as  the  point  for  which  the  distance  of  each  face  from  0  is 
proportional  to  the  surface  tension  of  that  face. 

Return  to  equation  2.1.44  to  determine  the  second  set 
of  equilibrium  conditions.  Differentiating  equation  2.1.35 
yields : 


(2.1.52) 


Substituting  2.1.52  into  equation  2.1.44  gives, 


(2.1.53) 


■ 
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Since 


/* 


and 


d  N  cannot  equal  zero. 


c*. 


(2.1.54) 


Thus  the  chemical  potentials  of  all  three  phases  are  equal. 
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2.2  Stability  Analysis 

For  the  purpose  of  stability  analysis  a  system  of  r 
components  with  a  spherical  nucleus  will  be  considered  to 
make  the  analysis  simpler.  Consider  the  two  systems  illus¬ 
trated  in  Figures  2.6(a)  and  2.6(b).  The  system  in  state 
A  is  a  large  volume  of  liquid  solution  at  pressure,  p\ 
surrounded  by  a  heat  bath  at  temperature,  T.  The  system 
shown  in  figure  2.6(b)  is  a  spherical  solid  nucleus  which 

is  close  to  its'  equilibrium  state.  The  nucleus  has  a  hydro- 

PM 

static  pressure  r q  .  The  radius  R,  of  the  nucleus  is  close 
to  the  equilibrium  radius, Rc .  The  total  number  of  moles 

in  the  nucleus, f/  ,is  close  to  the  number  of  moles  in  the 

a/ 

critical  nucleus. 

The  difference  in  free  energy  between  the  systems  is 


where 


and 

In  state  A 


(2.2.1) 

(2.2.2) 

(2.2.3) 


(2.2.4) 
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Constant 

Temperature 

Bath 


Figure  2.6(a) 


Figure 


Figure  2.6  Formation  of  a  Solid  Nucleus  in  a  Liquid  Volume 
where 

( 1 )  denotes  liquid  phase 
( 1 1  )  denotes  solid  phase 
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For  state  B  sufficiently  close  to  equilibrium^,  the  Helmholtz 
free  energy  can  be  expressed  as 

F6  =  ~  Ps'Vs'  '  PS  + 

r t^us.  /V  3  )  +  XD 

(2.2.5) 

The  system  in  states  A  and  B  is  constrained  by  the  fol¬ 
lowing  two  relationships: 


and 


(2.2.6) 


l  =  i}  i...  r, 

(2.2.7) 


The  pressure  in  state  A,  ,  is  not  significantly  different 

than  that  in  state  B,  ,  because  the  nucleus  is  small  rela- 

n> 

tive  to  the  volume  of  the  surrounding  liquid.  Since  the  pres¬ 
sure,  concentrations  and  temperature  of  the  bulk  liquid  phases 
in  states  A  and  B  are  constant,  the  chemical  potentials  are 
also  constant: 


hi 


. '  CP>  t 


,) 


/-/,  2...r.  (2.2.8) 


Assume  that  each  small  mass  element  is  in  a  state  of  local 
equilibrium  for  which  properties  are  the  same  function  of 
the  local  macroscopic  variables  as  at  the  equilibrium  state. 
"This  assumption  of  local  equilibrium  is  not  in  contradiction 
with  the  fact  that  the  system  as  a  whole  is  out  of  equilibrium" 
as  long  as  "the  dissipative  processes  are  sufficiently  domin¬ 
ant  to  exclude  large  deviations  from  statistical  equilibrium." 
(Glandorff  and  Progogine,  1971) 
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It  is  also  assumed  that  the  number  of  moles  in  the  inter¬ 
facial  zone  is  much  smaller  than  the  number  of  moles  in  either 
the  bulk  liquid  phase  or  the  solid  nucleus.  Thus  the  follow¬ 
ing  three  relationships  are  evident: 


and 


(2.2.9) 

(2.2.10) 

(2.2.11) 


Substituting  equations  2.2.4  and  2.2.5  into  2.2.1  yields 

A  F  -  -  (  f  8  '  ffl '  )  V £  -  (  ?8  -  P#  )  \Jq  '  +7?$  + 


r.  (2.2.12) 

Equations  2.2.9,  2.2.10,  and  2.2.11  are  used  to  simplify 
equation  2.2.12  to 

//  s 

AF  =  -aPVb  t ■  A~yc- A/  +&A  i  = 

(2.2.13) 

'/  I 

where  A  Pg  ~  y 


(2.2.14) 

//  / 

and  k*'*/ i  -  'Hjrt  ~  l  L  Z  2' ' '  C  , 

(2.2.15) 


Both  a/3 and  &  jy  (/- ',Z.  ,,c)  are  developed  for  the  near 


equilibrium  state. 
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Expanding  jP  and  a about  the  equilibrium  state  using 
a  Taylor  series  yields: 

AP(Z,n’)=&P(Pc,  ik-Plc-t  f  iA.P  .  AWS;  +...  ,  C- 1  ?.  ,.r , 

(2.2.16) 

(2.2.17) 

(2.2.18) 


where 


and 


AR-  R-R c 


4 i  = 


and 


A  A/ a 


V  <- 


t  &  R  t  db-S-j/' 

c)  R  Rcj^Cjt  ^  ^ j 

State  B  is  considered  to  be  sufficiently  close  to  the  equil- 


(2.2.19) 


R.Jltl  ^-r 


ibrium  size  nucleus  that  only  zero  order  terms  in  equations 

2.2.16  and  2.2.19  need  be  retained. 

1 

At  equilibrium  Laplaces  capillary  equation  (Defay  and  Pngo- 
gine,  1966)  is  used  to  relate  the  pressure,  radius  and  surface 
tension  of  a  spherical  nucleus.  From  equation  2.1.54. 


A  ^  (  Yd  j  NCjC  )  =  o  l  =/,2. . .  r 


(2.2.20) 


For  a  spherical  nucleus  the  volume  and  surface  area 
can  be  expressed  in  terms  of  the  radius  R,  which  is  near  the 


critical  radius,  R^, ; 


and 


V-  h/3  iY  R-2 


(2.2.21) 


1 


P u  -  P  ’  =  2  Y/R f 


A--  H fC 


(2.2.22) 
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Substituting  Laplace's  equation  2.1  and  equations  2.2.20 
to  2.2.22  into  equation  2.2.13  yields 

(2.2.23) 

A  plot  of  this  function  is  shown  in  figure  2.7  At  R=R 
the  free  energy;AF  ,  is  an  extremum.  At  the  critical  size  the 
conditions  of  equilibrium  are  satisfied.  It  is  obvious  from 
the  plot  of  this  function  that  the  critical  nucleus  is 
unstable.  A  nucleus  whose  radius  is  smaller  than  R^  will  com 
pletely  dissolve  and  a  nucleus  larger  than  the  critical  size 
will  spontaneously  grow.  This  theoretical  result  has  been 
confirmed  for  the  liquid-gas  bubble  equilibrium  (Tucker  and 
Ward,  1975).  The  implication  of  this  is  that  there  is  an 
energy  barrier,  equal  to  AF ,  to  the  nucleation  process 
This  is  the  energy  required  to  form  a  nucleus  of  at  least 
critical  size.  For  a  spherical  nucleus,  R  =  R q 

a  5*,,*  -  N/z  V'tt'  • 

(2.2.24) 

This  expression  can  be  extended  to  form  an  equation 
for  for  a  multifaceted  nucleus  in  the  following  way. 

Substitute  equation  2.2.20  into  equation  2.2.13  to  yield 

a  F  -  -a  P\i'a  +  r  A  • 


(2.2.25) 
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Figure  2.7  The  Variation  of  the  Change  In  Helmholtz  Free 

Energy  with  Radius 
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If  the  crystal  is  assumed  to  consist  of  planar  surfaces, 
(see  fig.  2.8),  V A  can  be  expressed  as  the  sum  of  surfaces 
of  the  nucleus: 


YAS= 

oC 

Similarly  it  is  possible  to  divide 
mids.  The  total  volume  is  then 


(2.2.26) 


the  volume  V  into  pyra- 


V*  |C(/3) 

The  change  in  pressure , A P  ,  is 


(2.2.27) 


A  P-  ( ZLY  )/h<K  (2.2.28) 

for  each  crystal  face  <=*»  .  The  change  in  free  energy,  obtained 
by  substituting  equations  2.2.26  to  2.2.28  into  equation 
2.2.25  is  equal  to: 

^  ^TOAx'  "  A  k  •  £ X  /O  (2.2.29) 

«< 

Equation  2.2.29  can  be  further  simplified  to 

°  rrnAx  ^  (2.2.30) 

This  relationship  will  be  important  in  section  3,1  concerning 


nucleation  theory. 
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Figure  2.8  Theoretical  Nucleus  Shape  Consisting  of  Planar 
Surfaces 
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2.3  Expression  for  Equilibrium  Size  of  Nucleus 

It  is  now  important  to  derive  an  expression  for  the 
equilibrium  size  of  the  nucleus  in  terms  of  the  measurable 
quantities  temperature  and  pressure.  The  system  described 
in  section  2.2  will  be  used  with  the  additional  restriction 
that  it  consist  entirely  of  one  component.  At  equilibrium 


4 


C  T,  T,  P") 


(2.3.1) 


^  "CTjP")  "C  P')  . 

Substituting  equation  2.3.1  into  equation  2.3.2 

^  ’C  "(  t>  p‘>  -  p"-p) 

T0  ,  the  bulk  equilibrium  freezing  temperature, 
ature  for  which 


(2.3.2) 

yields 


(2.3.3) 

is  that  temper- 


//  / 

st/C  /a>  (7~a)p'>.  (2.3.4) 
The  chemical  potentials  in  equation  2.3.3  can  be  expressed 
in  terms  of  T  by  using  the  Gibbs-Duhem  relation: 


T~ 

^Ct,  p)--~Y(Te,p->  -J  sacr 


(2.3.5) 


Thus 


sy'C  )  ;  P')  -  P  )  ~ To ^  ^  d-T  (2.3. 

P')-fs'eL 7(2>3> 


6) 


7) 


and 
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Substituting  equations  2.3.6  and  2.3.7  into  equation 
2.3.3  and  using  the  equality  2.3.4  yields: 

-S is'-s'^cS-r^v-"  . 

To  (2.3.8) 

The  latent  heat  of  fusion, L,  is  defined  as, 

L  =  T(s'-s")  . 

(2.3.9) 

Substituting  equation  2.3.9'  into  equation  2.3.8  yields 

Jr*  lt-Cp'-p') 

(2.3.10) 

Let  <  be  the  average  latent  heat  over  the  temperature 

range  T^to  T  and  assume  that(L/T)  is  not  strongly  dependent 
on  temperature.  Integrating  equation  2.3.10  over  the  temper¬ 
ature  range  T^to  T  yields: 

Cf'-P')-  <u~?  In(r°/r)  A" . 

(2.3.11) 

Substituting  this  value  of  A  p  into  equation  2.2.28  yields: 

(V*=  ay  (%rT)).  <2.3.12) 

where  h  is  the  equilibrium  distance  of  face  from  point  O. 
Note  that  h  is  a  function  of  temperature  and  pressure,  the 
pressure  dependence  entering  via  the  dependence  of  To  on 
pressure.  This  is  evident  in  Figure  2.9  which  is  a  phase 


diagram  of  a  substance  that  expands  on  solidification.  As 
the  pressure  increases,  the  temperature  at  which  the  phase 


Figure  2.9  Phase  Diagram  for  a  Typical  Liquid  that  Expands 
on  Freezing. 
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transition  occurs  increases  as  well.  L  depends  on  pressure 
to  a  lesser  extent. 

Figure  2.10  is  a  phase  diagram  for  a  typical  cooling 
process  at  constant  pressure.  Liquid  at  pressure  P  can 
be  cooled  to  point  A,  the  bulk  freezing  point,  and  all  points 
along  this  curve  are  in  stable  equilibrium  states.  The  por¬ 
tion  of  the  curve  from  A  to  B,  the  spinodal  point,  is  meta¬ 
stable.  That  is,  the  system  will  remain  in  the  original 
state  as  long  as  no  large  property  perturbations  occur. 
Microscopically  the  substance  is  always  fluctuating  in  density 
and  if  a  disturbance  large  enough  to  take  a  portion  of  the 
system  to  state  D  occurs  that  part  of  the  system  will  remain 
as  state  D.  This  implies  that  there  is  a  critical  size  of 
a  nucleus.  The  smallest  amount  of  material  that  remains 
in  state  D  is  called  the  critical  nucleus.  Under  a  specified 
temperature  and  pressure  it  is  possible  with  nucleation  theory 
to  estimate  the  frequency  of  nucleus  formation  in  a  liquid. 

The  next  chapter  will  treat  this  subject  in  more  detail. 


Temperature 
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T0 
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Figure  2.10  Phase  Diagram  Illustrating  Cooling  Process  at 
Constant  Pressure 


3.0  Theory  and  Calculations 

The  theoretical  developments  derived  in  chapter  two 
from  equilibrium  thermodynamics  are  applicable  to  any 
liquid-solid  equilibrium.  In  this  chapter  a  pure  ice-water 
system  is  used  to  illustrate  how  these  relationships  can  be 
used  to  determine  the  anistropic  surface  tensions  for  ice. 

A  nucleus  shape  related  to  the  shape  of  the  unit  cell  is 
assumed.  Figure  3.1,  is  a  composite  diagram  of  experimental 
observations  of  the  shape  of  the  unit  cell  at  varying  temper¬ 
atures  and  pressures.  From  this  figure,  it  can  be  seen  that 
for  ice  at  1  atm  and  temperature  between  0°C  and  -50°C  the 
unit  cell  is  a  hexagonal  prism,  commonly  known  as  Ice  I ^  . 

For  this  shape  two  distinct  surface  tensions  are  cal¬ 
culated,  y,  ,  corresponding  to  the  basal  surface  (perpen¬ 
dicular  to  the  c  axis)  and  2^  ,  corresponding  to  the  pris¬ 
matic  surfaces  (perpendicular  to  the  a  axis),  see  Figure  3.2. 
Using  nucleation  rate  theory  and  experimental  results  for  the 
maximum  supercooling  a  relationship  between  the  surface  ten¬ 
sions  of  the  two  different  surfaces  of  the  water  nucleus  is 
developed.  Following  a  short  discussion  of  contact  angle 
theory,  experimentally  obtained  data  for  the  surface  energy 
of  the  basal  plane  of  ice  are  used  to  determine  unique  values 
for  and  2^  . 
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Temperature  (°C ) 
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Figure  3.1  Phase  Diagram  of  the  Solid  Phases  of  Water 
(adapted  from  Hobbs,  1974) 
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Fiqure  3.2  Shape  of  unit  cell  for  Hexagonal  Ice 

Where  h(  =  perpendicular  distance  from  centroid 
to  basal  plane 

h  2  =  perpendicular  distance  from  centroid 
to  prismatic  plane. 
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3 . 1  Nucleation  Rate  Theory 

Small  water  droplets  can  be  cooled  well  below  0°C  without 
initiating  freezing.  This  is  because  small  droplets  are 
less  likely  to  be  contaminated  with  foreign  particles  to 
catalyze  heterogeneous  nucleation.  Homogeneous  or  spontaneous 
nucleation  will  occur  as  the  liquid  is  cooled  below  the  bulk 
freezing  temperature  (0°  for  water.)  In  supercooled  water 
there  are  microscopic  fluctuations  in  density, temperature 
and  pressure  ( Frenkle, 194 6 ) .  Thus  tiny  ice  nuclei  are  con¬ 
stantly  being  formed.  In  chapter  two  it  was  shown  that  if 
a  nucleus  of  a  "critical"  minimum  size  is  created  it  will 
spontaneously  grow.  Since  it  is  impossible  to  measure  the 
probability  of  a  nucleus  being  formed  in  the  supercooled 
liquid,  homogeneous  nucleation  theory  was  developed.  This 
theory  predicts  the  number  of  critical  nuclei,  J,  which  are 
spontaneously  formed  in  a  unit  volume  of  liquid  per  unit 
time.  Becker  (1938)  investigated  liquid-solid  systems  and 
proposed  the  following  form  for  the  equation  for  the  rate 
at  which  nuclei  are  formed: 

T-  ~A  F J(k  t)  )  ; 

(3.1.1) 

where  Z  is  the  rate  constant  and  K  is  Boltzmann's  constant. 
Turnbull  and  Fisher  (1949)  expanded  this  theory  and  developed 
the  following  expression  for  the  rate  constant ; ^  , 
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2:=n  KT/h  , 

7  1  (3.1.2) 

where  n  is  the  number  of  nuclei  per  unit  volume  and  h  is 
Planck's  constant. 

Turnbull  (1952)  generated  dispersions  of  small  droplets 
and  used  a  dilatometer  to  experimentally  demonstrate  the 
difference  between  heterogeneous  nucleation  and  homogeneous 
nucleation.  Values  for  surface  tension  calculated  by  his 
technique  are  for  spherical  nuclei.  Hesstvedt  (1964)  supercooled 
small  water  droplets  and  using  homogeneous  nucleation  theory 
calculated  a  value  of  31.7  dyne/cm  for  the  interfacial  surface 
tension  at  0°C.  Hesstvedt  (1964)  derived  the  following  equa¬ 
tion  to  calculate  Z: 

2S  v/3  irr3n  xt.  0.1.3) 

where  r  is  the  radius  of  the  water  droplet.  From  Figure  33 
it  is  seen  that  a  one  micron  water  droplet  will  spontan¬ 
eously  freeze  at  -38.8°C.  This  experimentally  determined 
maximum  supercooling  is  assumed  to  be  the  homogenous 

nucleation  temperature.1 *  Using  these  values  and  substituting 

24 

into  equation  3.1.3,  Z,  the  rate  constant  is  6.27  x  10 
nuclei/cm  /s.  With  this  value  of  Z,  and  assuming  a  value 
of  J=l,  AF  can  be  expressed  in  terms  of  the  critical 

f  WA  A 

1This  value  is  very  close  to  the  homogeneous  nucleation  tem¬ 

perature  reported  by  other  researchers  eg:  Schaefer  (1948) 

reports  -38.9°C,  Schaefer  (1952)  reports  -38.5°C,  Weickman 
(1949)  reports  (-40°C). 


Temperoture  (°C) 
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Figure  3.3  Relationship  Between  Water  Droplet  Size  and 
Homogeneous  Nucleation  Temperature. 

(Adapted  from  Hesstvedt,  1964) 


46 


size  based  on  the  results  in  the  last  chapter.  This  equation 
can  then  be  used  to  predict  the  supercooling  if  the  surface 
tensions  are  known  or  (as  illustrated  in  the  following  section) 
to  predict  the  relationship  between  the  surface  tensions 
given  the  supercooling. 
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SL.  SL 

3.2  Calculation  of  Relationship  of  ~6 1  and 

Si-  SL 

To  calculate  the  relationship  between  *2f,  and  for  the 
ice-water  system  it  is  necessary  to  assume  a  nucleus  shape. 

Most  researchers,  like  Turnbull  (1950),  assume  that  the  nucleus 
is  spherical  in  shape.  Hesstvedt  (1964)  realized  that  unique 
faces  of  a  crystal  should  have  unique  surface  tensions,  and 
that  the  nucleus  is  not  necessarily  spherical  but  he  calcu¬ 
lated  a  weighted  mean  value  for  different  shapes.  As  discussed 
in  the  introduction  to  chapter  three  in  this  work,  it  will 
be  assumed  that  the  nucleus  shape  is  related  to  the  structure 
of  the  unit  cell  at  the  nucleation  temperature  and  pressure. 

It  is  clear  that  the  surface  area  of  the  crystal  can  be  for¬ 
mulated  in  terms  of  h  j  and  h^  see  Figure  3.2,  the  perpendi¬ 
cular  distance  from  the  centre  of  the  crystal  to  the  basal 
and  prismatic  planes , respectively .  The  ratios  of  h  ,  and  hj. 

not  be  the  same  as  in  the  unit  cell.  Henceforward  quanti¬ 
ties  subscripted  as  1  refer  to  the  basal  plane  and  those 
subscripted  as  2  refer  to  the  prismatic  plane. 

The  total  area  of  the  two  basal  planes,  A  (  ,  is; 

CofckO  ,  (3.2.1) 

The  total  area  of  the  six  prismatic  faces,  ,  is; 

Ax  -  C-O+  (c0° 


(3.2.2) 
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It  is  possible  to  express  the  distances  h|  and  h^_in 
terms  of  the  surface  tensions  of  their  respective  planes, 
2^  SU  and  ^ 


and 


(as  seen  in  chapter  2); 

h,-ZLY,  c" I  (<  L“7\n  C  ) , 

In ,2  <s" / [ ^L'?ln(7b^) 


(3.2.3) 


(3.2.4) 


By  substituting  equations  3.2.3  and  3.2.4  for  h  j  and 
h  ^  into  the  area  formulations,  expressions  relating  the  area 
to  surface  tension  and  temperature  are  obtained; 

')  J 


jun 


L 


In  Vfolr) 


and 


y  In  Cr°/f')  - 


3.  SL  St- 

.  co-tO’O0 .  ^ 


(3.2.5) 


(3.2.6) 


Substituting  equation  3.2.5  and  3.2.6  into  equation  2.2.30 
yields;  7-  a.  e  -r  /  \ 

2  Cr'/  CofleO  /(<L>)n(  t))  ■ 


(3.2.7) 


Equation  3.2.7  can  be  substituted  into  3.1.1  the  nucleation 
rate  equation  to  give; 

J  =  2  e*P  (-■<}-?  ir"  co’t(eO°J((<L>lr1(r°/1))  kr)  . 

(3.2.8) 


The  exponential  nature  of  equation  3.2.8  suggests  that  the 

Sl 

relationship  between  o,  and  is  relatively  insensitive  to 
the  values  of  J  and  Z.  For  the  purpose  of  calculation  J, 
the  nucleation  rate,  is  set  equal  to  1.  The  rate  constant 


49 


used  will  be  the  one  derived  from  Hesstvedt  (1964)  in  section 


S 54 

3.1.  Solving  equation  3.2.8  for  2j  )using  the  properties 

of  water  listed  in  Appendix  I  ,  yields; 


(3.2.9) 


Figure  3.4  is  a  graph  of  this  relationship.  A  graph 

-v  S4 

of  the  percentage  difference  in  o-l.  for  two  assumed  values 
of  J  (  J=1 ,  and  J=100)  is  presented  in  Figure  3.5.  .  A  differ¬ 
ence  in  the  assumed  rate  of  two  orders  of  magnitude  of  10 

causes  only  a  4%  difference  in  .  Thus  it  is  seen  that 

SL  . 

Ya.  12  relatively  insensitive  to  the  assumed  value  of  J. 


Si- 


In  the  following  section  unique  values  for  7t,  and  Ox. 
are  determined  by  using  contact  angle  theory  and  experimental 

data  to  arrive  at  a  value  for  the  surface  tension  of  the 

.  SL 

basal  plane  of  the  ice-liquid  water  interface, 
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Y  SL 

0 1  (dyne /cm  ) 


Relationship  Between  the  Surface  Tensions  of 
the  Basal;  JS^and  Prismatic^TT^  5  Planes  of  Ice 


Figure  3 . 4 
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Figure  3.5  Effect  of  Varying  the  Assumed  Rate  at  Which  Nuclei 
are  Formed  on  . 
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3.3  Contact  Angle  Theory 

SL  SL. 

In  order  to  determine  unique  values  for  and  "7fx  , 
contact  angle  data  for  carbon  disulphide  (CS^  )  droplets  on 

the  basal  plane  of  ice  obtained  by  Adamson  (1970)  are  used. 

.z.s/ 

Equation  of  state  relationships  between  o  ,  liquid-vapour 

v  SV 

0  ,  solid-vapour,  and  o  ,  solid-liquid  surface  tension 

(Ward  and  Neumann,  1974)  derived  by  Ward  are  used  to  correl¬ 
ate  Adamson's  data  to  the  basal  plane  of  the  ice-liquid  water 
system . 

The  angle  at  which  a  liquid  wets  a  solid  (Figure  3.6) 
when  the  three  phases  are  in  equilibrium  is  ,  the  contact 
angle.  The  contact  angle  is  related  to  the  surface  tensions 

between  the  three  phases  by  Young's  (1805)  equation; 

cv  SL 

(3.3.1) 

Ward,  (1975)  states  that  "in  any  system  in  which  a  smooth, 
homogeneous  and  inert  solid  is  in  equilibrium  with  a  liquid 
and  a  vapour  phase,  adsorption  must  take  place  at  the  solid- 
vapour  interface  and  the  adsorption  must  be  positive." 

Whalen  and  Wade  (1967)  and  Whalen  (1968)  measured  the  adsorp¬ 
tion  of  liquid  octaine  on  Teflon  (contact  angle  26°)  and 
found  that  the  solid-vapour  adsorption  was  significant. 
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Figure  3.6  Three  Phases  in  Equilibrium  with  Each  Other. 

(Adapted  from  Defay  and  Progogine,  1966) 
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Using  the  system  shown  in  Figure  3.7  and  the  follow¬ 
ing  assumptions; 

1.  solid  is  non-dissolving 

2.  solid  has  negligible  vapour  pressure  so  that  there 
is  a  single  component  present  in  the  liquid  and 
vapour  phases. 

3.  solid  is  smooth,  homogeneous  and  rigid, 

Ward  ensures  that  a*  is  the  equilibrium  contact  angle.  Using 
the  method  of  Gibbs  (1961)  for  equilibrium  of  this  system, 


1.  uniformity  of  temperatures 

_  3l  sk/ 

=  -Y  =-y 

sl  sv  lv/ 

Ward  shows  that  ^  is  a  function  of  if  and  V  ie; 


and  2 . 


SL  .  5V  lV\ 

«  FC-u  ,  v  ) 


(3.3.2) 


To  develop  the  explicit  relationship  of  equation  3.3.2  the 
following  derivation  follows  that  of  Ward  (1975).  To  deter¬ 
mine  this  function  the  following  variables  will  be  used: 


nl 


-  Y 


(3.3.3) 


and 


L.V  S'/ 

£L-  -6  +  i 


(3.3.4) 

4.1/  SV 

and  can  be  expressed  in  terms  of  the  variables 


and 


ie ; 


and 


L.V 


(/7£  +  £  ) 


(3.3.5) 


-i) 


(3.3.6) 
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Figure  3.7  Three  Phase,  Two  Component  System 
(Adapted  from  Ward,  1975  ) 


iiOTtBBBawaBaBBg 
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Substituting  equations  3.3.5  and  3.3.6  and  3.3.1  into 
equation  3.3.2  yields: 

E')  ~  F (&,  /p  )  -  ((^  ^i)oos 

(3*3*8) 
^4 


where  7  is  a  function  of  E  and 


/rt  because  it  is  a  function 
of  o  and  o  .  These  equations  are  subject  to  the  follow¬ 


ing  three  limits: 

$v 

1.  if  ^  ~ 

then  COSO  must  equal  one  and  therefore 

F(f  o)  =o 

SV 


(3.3.9) 


(3.3.10) 


2.  as  if  approaches  zero,  ^  approaches  TT*  ,  therefore 


F(/f,  /f  )-  ST? 


(3.3.11) 

LV  S’U 

3.  as  ^  approaches  zero,  and  remains  finite,  the 
contact  angle  will  approach  zero  and  therefore 


F j  )- 


(3.3.12) 


When  adsorption  is  present: 


and  thus 


z.1/ 

y  ^  y 

-  C? 


(3.3.13) 


SL 


(3.3.14) 

Therefore,  ,  a  function  of  ^  and/^  can  be  expanded  in  a 
Taylor  series  about  F  {^o)  as  follows; 


2 (3.3.15) 


sVl-O 
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where  higher  order  terms  are  neglected. 
After  defining  the  following  functions. 


and 


'■2  A  rvf  ‘fa  - 


(3.3.16) 


(3.3.17) 


<)  (V(J-  r  O  r 

and  using  equation  3.3.10,  equation  3.3.15  can  be  written  as: 


(3.3.18) 

Since  equation  3.3.18  must  also  satisfy  the  limits  of  3.3.11 
and  3.3.12, 


and 


3(t)  =  0 

■f  (i)-rVi 


(3.3.19) 


(3.3.20) 


Substituting  back  into  equation  3.3.2  yields  the  following 
equation : 


SL 

X- 


(3.3.21) 


Equation  3.3.21  is  the  equation  of  state  in  the  presence 
of  adsorption.  Equation  3.3.21  can  be  combined  with  equation 
3.3.1  to  yield  another  relationship  which  will  be  useful 
in  the-  following  analysis. 


S'*  LV  ^ 

ar  -  ( 


(3.3.22) 
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In  the  following  paragraphs,  a  value  of  (basal  plane 
of  ice  -  CS^  vapour)  is  obtained  using  the  previously  devel¬ 
oped  contact  angle  equations  and  Adamson's  data.  This  value 
SV  SV 

of  is  assumed  to  be  approximately  equal  to  ^  (basal 

plane  of  ice-water  vapour). 

Experimental  evidence  indicates  that  the  solid-vapour 
interfacial  tension  is  not  strongly  dependent  on  the  adsorp¬ 
tion  of  the  contacting  vapour  or  the  liquid-vapour  surface 

tension.  From  contact  angle  data  Hellwig  and  Neumann  (1968) 

SV 

found  that  V  the  solid-vapour  surface  tension  is  nearly  con¬ 
stant  for  various  fluids  on  a  solid. 


s  l 

The  equation  of  state  is  then  used  to  calculate  # 

SL 

(basal  plane  of  ice-liquid  water).  Knowing  K  for  the  basal 

SL 

planes  of  ice, ^  for  the  prismatic  planes  can  be  calculated. 

Figure  3.8  is  a  graph  of  the  contact  angles  that  Adamson 
(1970  )  measured  for  different  temperatures  of  liquid  carbon 
disulphide  on  the  basal  plane  of  ice.  At  -38.8°C  the  contact 
angle  for  ice  is  about  29°C.  A  linear  regression  for  CS  2. 
liquid  vapour  surface  tension  yielded  the  following  equation: 


LV  I'lJL. 

-  7-f.  (/.-  'r/s'4-9 .)  (3.3.23 

Sv 

Using  Adamson's  data,  and  equations  3.3.22  and  3.3.23,  $  ; 

the  surface  tension  for  the  basal  plane  of  ice  -  CS^  vapour 

is  calculated.  At  the  measured  supercooling  of  -38.8°C,  the 

.sv 

calculated  value  for  if  is  36.1  dyne/cm. 
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Figure  3.8  Contact  Angle  for  Carbon  Disulfide  on 
Basal  Plane  of  Ice. 

(Adapted  from  Adamson , 19 7 0 ) 

where  "0"  cooling  sequence 
"t"  heating  sequence 
"8"  cooling  sequence  (run  2) 

"0"  heating  sequence  (run  2) 
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sv 


This  value  is  assumed  to  be  equal  to  (basal  plane 

of  ice-water  vapour). 

sv 

To  use  equation  3.3.21  to  calculate  If  (basal  plane 

4  V 

of  ice-water)  a  value  of  If  (liquid-water  vapour)  must  be 
determined.  This  was  done  by  using  Dorsey's  (1968)  data 

J  ^  ^  f  r4  T 

for  the  variation  of  &  V  J  ^  with  temperature; 


d'i'/ct'T  -  qT  +-  b 


(3.3.24) 


where  T  is  C. 


From  Dorsey  (1968)  the  value  for  a  is  -0.006  dyne/cm  -  °C/ 


C  and  the  value  for  b  is  -0.139  dyne/cm  -  °C. 

Integrating  equation  3.3.24  between  20°C  and  -38.8°C  yields: 


-3%-  $ 


r 


d  V  -  (  ,  c>  o  e>  Lj>  T  -+  c> .  i  3  ?  ^  c| 


2D. 


J 
20. 


(3.3.25) 


This  equation  can  be  further  simplified  to; 

-3  8  ?  -  58. 8 

If  ~  V  ~  ,  o  oo3  T  j  -  o  ./ 39T  I 

-3?.?  10.  zo  10  (3.3.26) 

Using  the  value  of  12. If  dyne/cm  (from  Dorsey,  1968)  at 

20°C,  the  liquid  vapour  surface  tension  at  -38.8°C  is  80.6 

dyne/cm.  This  value  agrees  with  the  value  obtained  from  the 

equation  in  the  A.S.M.E.  steam  tables  (1977). 

SV  LV 

Substituting  the  values  of  If  and  ~6  into  equation 
3.3.21  gives  a  vlue  of  17.0  dyne/cm  for  the  basal  plane  of 
ice-water  interfacial  tension.  By  using  equation  3.3.9  from 
the  previous  section  a  value  of  28.9  dyne/cm  is  calculated 


61 


for  the  liquid  solid  interfacial  tension  for  the  prismatic 
planes.  The  values  of  h  j  and  h  2,  calculated  from  equations 
3.2.3  and  3.2.4  are  shown  in  Figure  3.9. 

As  shown  earlier  the  values  of  surface  tension  for  the 
balas  and  prismatic  planes  are  insensitive  to  the  chosen 
value  of  J,  the  rate  at  which  the  nuclei  form.  It  is  also 
evident  that  these  values  will  be  equally  insensitive  to 

changes  in  Z,  the  rate  constant.  A  6%  change  in  the  calcu- 

ri  sl  0 

lated  value  of  '?fl  will  change  ^  by  3%.  A  1C  change  in 

o  SL 

the  estimated  subcooling  of  -38.8  C  will  change  K,  by  1%  and 
SL  Si-  SL 

by  3%.  The  values  of  z,  and  are  very  sensitive  to 

L * 


the  estimated  liquid-vapour  surface  tension.  A  change  in  K 

s  u 

of  only  5%  causes  a  36%  change  in  the  calculated  value  of 

S  L 

and  a  25%  change  in  the  calculated  value  of 


62 


Figure  3.9  Diagram  of  Proposed  Ice  Crystal  Nucleus 
(drawn  to  relative  scale) 


4 . 0  Discussion 


Evidence  by  Simpson  et  al.  (1974a,  1974b)  and  Hillig 
(1958)  indicates  that  there  is  an  orientational  dependence 
of  the  surface  tension  for  ice-water  systems.  Kallungal 
(1975)  believes  that  "since  ice  crystals  show  definite  pre¬ 
ferred  orientation  in  growth  the  it  associated  with  each  of 
the  growth  directions  should  be  considerably  different." 

The  values  of  17.0  dyne/cm  for  the  basal  plane  and  28.9 
dyne/cm  for  the  prismatic  plane  (both  at  -38.8°C)  are  consist¬ 
ent  with  what  Ketchum  and  Hobbs  (1969)  observed.  They  noted 
that  "when  the  basal  plane  was  exposed  at  a  grain  boundary, 
the  groove  was  generally  asymmetrical,  and  the  results  indi¬ 
cated  that  the  surface  energy  of  the  basal  plane  was  slightly 
lower  than  that  of  the  other  crystallographic  faces"  (Hobbs, 
1974).  Higuchi  (1961)  also  found  that  "the  basal  plane  prob¬ 
ably  has  the  lowest  value  of  the  surface  energy  of  all  planes 
in  ice"  (Hobbs,  1974).  Ketchum  and  Hobbs  (1969)  note  that 
only  3%  of  their  total  number  of  datum  points  include  cases 
where  the  basal  plane  was  exposed  at  the  grain  boundary  groove. 
This  means  that  their  value  of  (29-3)  dyne/cm  at  0°C  which 
they  extrapolate  to  (25  -  3 ) dyie/cm  at  -40°C  is  primarily  that  of 
the  prismatic  plane.  This  value  agrees  well  with  the  value 
of  28.9  dyne/cm  obtained  in  this  work. 

Experiments  done  by  Ketchum  and  Hobbs  (1967)  indicate 
that  an  ice  crystal  will  tend  to  grow  in  a  direction 
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perpendicular  to  its  c-axis  in  subcooled  water.  These 
results  are  corroborated  by  experiments  conducted  by  Fernandez 
(1967  a)  Poisot  (1968),  Vlahakis  (1972),  Simpson  et  al.  (1974 
b)  and  Huige  and  Thyseen  (1969).  Ketchum  and  Hobbs  (1967) 
propose  that  a  direction  parallel  to  the  basal  plane  is  the 
preferred  direction  of  growth  for  a  crystal  because  it  will 
maximize  the  low-energy  basal  surface.  This  is  consistent 
with  the  values  obtained  in  this  paper. 

A  disadvantage  of  this  method  is  that  it  yields  values 
of  "if,  and  at  only  one  temperature  -  that  of  the  homogen¬ 

eous  nucleation  point.  To  compare  the  values  at  -38.8°C  to 
those  of  other  researchers  at  0°C  it  is  necessary  to  determine 


s  u 

how  *2^,  and  2fz  vary  with  respect  to  temperature  ie: 


(4.0.1) 


and 


(4.0.2) 


in  order  to  extrapolate  the  surface  energies. 

Figure  4.1'*'  is  a  graph  of  the  variation  of  Jd~T  with 


temperature.  From  this  graph  it  is  seen  that, 

J  2/  /  9  ~J~  rs  0.132  dyne/cm/°C  (4.0. 

SL 

Isolating  "2^  from  equation  3.2.8  and  differentiating  with 
respect  to  temperature  yields 


^The  points  on  the  graph  are  calculated  from  Adamson's 
contact  angle  data  at  various  temperatures,  Ward's  equations 
of  state  and  Young's  equation. 


(  dyne /cm  ) 
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Figure  4.1  Variation  of  Liquid-Solid  Interfacial  Surface 

Tension  for  the  Basal  Plane  of  Ice  with  Temperature 
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}  *1/0  7  =  <p  [tin  (  VT;  -  pQ  /^S- 

~  £  Int^/r) /j’J  /l j  >  (4.o.4) 

where 

(p  -  In  (?:)  K  ^L.~p  /(/{.  %  cc>  ■[  (pC?  if"  ) )  (4.4.5) 

and  the  temperature  dependence  of  In  ^2:,  K,  <d  L?,  and  v  ^  is 
ignored.  At  38.8°C 


SL 


dl*  /JT  =  -.18 

(4.0.6) 

Thus 

o  ,  J 

at  0  C  the  values  of  ot  and 

.SL 

are : 

Y  SC 

V  =  22. 

0  dyne/cm 

(4.0.7) 

and 

y  54 

=0.0 

dyne/cm 

(4.0.8) 

(no  neg  values  allowed). 

.SL 

These  values  for  "3,  the  solid-liquid  surface  tension,  are 
not  intended  to  be  presented  as  accurate,  because  the  calcu¬ 
lated  values  for  the  slopes  are  valid  at  -38.8°C  but  are  not 
necessarily  at  0  C.  The  slope  of  c>^l<s~T  j_s  larger  than 
Hobbs'  (1974:441)  suggested  range  of  0.1  to  0.35  dyne/cm/°C. 

SL.  ___ 

It  is  however  interesting  to  note  that  /&  >  is  positive 

$c  .  SL 

and  ^  J  ^  is  negative.  This  indicates  that  the  o, 

,$L  S/~ 

values  of  o ,  and  “B-z.  could  approach  each  other  as  the 
temperature  approaches  0°C.  This  could  explain  why 
most  researchers  do  not  notice  any  orientational  dependence 
for  solid-liquid  surface  tensions  at  0°C.  It  is  suggested 


that  the  anisotropy  of  the  solid  liquid  interfacial  tension 
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becomes  more  prominent  as  temperature  decreases. 

Further  experimentation  is  required  to  verify  the 
results  proposed  in  this  paper.  In  general  there  is  an  orien¬ 
tational  dependence  on  surface  tension  in  ice-water 
systems  and  the  basal  plane  (as  expected)  has  the  lower  sur¬ 
face  tension. 

The  approach  here  is  equally  applicable  to  other  liquid- 
solid-vapour  systems. 
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1. 


2  . 


3  . 


4.1  Conclusions 

SL 

There  is  anisotropy  of  ^  in  the  ice-water-water  vapour 
system  at  -38.8  °C. 

Using  equilibrium  theory, and  homogeneous  nucleation 
theory  and  experimental  data, it  is  possible  to  generate 
a  relationship  between  the  solid-liquid  surface  tensions 
of  an  assumed  nucleus  shape. 

For  ice-water  this  relationship  is 

%  (  K.  )=  1.  A  x  104  (dyne3/cm3) 

(4.1.1) 

4.  For  ice-water  at  -38.8°C. 

31 

X  =17.0  dyne/cm 


S  4 

y2  = 


28.9  dyne/cm 


(4.1.2) 

(4.1.3) 
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Appendix  I:  Properties  of  Water 


L  =  333.475  x  107  ergs/g 

v  =  1.0906  cm2 /g 

T0  =  273. 2°K 

T  =  234.4°K 

K  =  1.38054  x  10~16  erg/°K 

h  =  6.62554  x  10”27  erg/s 

n  =  3.067075  x  1022  molecules 
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Appendix  II: 


Geometric  Relationship  of  Area  in  Crystal 


Euler's  relationship  between  the  area  and  the  linear 
dimensions  is  a  second  order  homogeneous  function;  thus, 


/)y3( A V), , . . .  \  1^0  =  A 

A .  II. 1 

Differentiating  this  equation  with  respect  to  A.  yields 

}  i)  cT\  JZX/i/v)  A.  1 1. 2 


Since  this  is 
ing  the  terms 


true  for  all 
yields : 


is  set  equal 


to  1.  Rearrang 


A.  II .  3 
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